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RESUME 

On  etudie,  a l'aide  de  solutions  numeriques,  l'etat  stationnaire 
de  la  defocalisation  thermique  d'un  faisceau  laser  a impulsions  multiples 
soumis  a un  vent  lateral.  On  trouve,  pour  un  taux  d' impulsions  favorable, 
une  augmentation  de  100%  de  l'intensite  maximale  par  rapport  au  faisceau 
qui  se  propagerait  sans  deformation  mais  avec  absorption  lineaire.  La 
valeur  correspondante  du  gain  exprime  en  terme  de  l'intensite  moyenne  sur 
la  surface  de  la  cible  qui  regoit  85%  de  la  puissance  totale  est  de  40%. 
Trois  parametres  de  similitude  suffisent  pour  definir  compl&tement  le 
Probleme:  le  nombre  k qui  est  egal  a la  reciproque  du  nombre  d' impulsions 

transmises  pendant  le  temps  de  transit  de  l'air  dans  le  faisceau,  le 
nombre  de  Fresnel  F base  sur  la  longueur  caracteristique  de  1 ' interaction 
absorption-intensite  et  le  nombre  y relie  a 1' absorption  lineaire.  Une 
etude  parametrique  en  fonction  de  ces  trois  parametres  demontre  que  des 
conditions  optimales  de  gain  dans  le  cas  d'un  faisceau  gaussien  collimate 
existent  pour:  0.8  s k < 1.0  et  F > 20.  L'amplification  decroit  gene- 

ralement  avec  1 'accroissement  de  y mais  de  fa<jon  peu  sensible.  (NC) 


ABSTRACT 


Steady  state  thermal  blooming  of  multipulse  laser  beams  sub- 
jected to  cross  winds  is  studied  in  the  light  of  numerical  solutions. 
Enhancement  resulting  in  a 100%  increase  of  the  peak  irradiance  with  re- 
spect to  the  non-bloomed  but  linearly  absorbed  beam  is  demonstrated  at 
the  optimal  pulsing  rate.  The  corresponding  gain  in  terms  of  spatially 
averaged  power  density  over  the  target  area  containing  85%  of  total 
power  is  40%.  A parametric  study  is  made  as  a function  of  three 
similarity  parameters  which  are  sufficient  to  completely  define  the 
problem:  the  number  jc)  which  is  the  reciprocal  of  the  number  of  pulses 

per  flow  time,  the  Fresnel  number  F based  on  the  absorption-irradiance 
interaction  length  scale  and  the  absorption  number  y related  to  linear 
depletion.  Optimal  conditions  for  a collimated  Gaussian  beam  are  shown 
to  exist  for  0.8  -</ 1.0  and  F >)  20.  Intensification  generally  de- 
creases with  increasing  y but  only  slightly.  (U) 
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1.0  INTRODUCTION 

Atmospheric  applications  of  high-power  laser  beams  are  subjected 
to  the  nonlinear  distorting  mechanism  called  thermal  blooming.  As  a result 
of  absorption  of  a fraction  of  the  beam  power,  the  enthalpy  of  the  air  is 
slightly  increased  along  the  propagation  path.  This,  in  turn,  alters  the 
distribution  of  the  refractive  index.  Although  these  changes  are  relatively 
small  owing  to  the  low  absorption  coefficient  of  the  atmosphere  at  the  wave- 
lengths of  interest,  they  can  nevertheless  give  rise  to  significant  distortion 
of  the  laser  beam  irradiance  profile  at  ranges  of  a few  kilometers. 

The  thermal  blooming  process  has  been  the  subject  of  extensive 
theoretical  studies  but,  because  of  the  complexity  of  the  modelling  differential 
equations,  computer  programs  had  to  be  developed  to  obtain  solutions.  Exper- 
iments were  mostly  performed  in  the  laboratory  where  the  thermal  blooming 
phenomenon  was  enhanced  in  a controlled  atmosphere.  A comprehensive  review  of 
the  work  done  on  this  subject  may  be  found  in  Refs.  1-5.  Investigations 
on  the  thermal  blooming  process  were  also  undertaken  at  the  Defence  Research 
Establishment  Valcartier  (DREV) . Refs.  6-8  present  a theoretical  model 
with  emphasis  on  the  definition  and  the  use  of  similarity  parameters. 

Ref.  9 describes  a simulation  experiment  and  Ref.  10  summarizes  the 
principal  features  of  a computer  program  based  on  the  DREV's  theoretical  model. 


This  report  describes  a study  of  the  thermal  blooming  of  multipulse 
laser  beams  using  computer  solutions.  A reduction,  or  even  an  elimination, 
of  the  adverse  thermal  blooming  effects  can  be  expected  if  high-power  lasers 
are  operated  in  a multipulse  configuration  rather  than  in  a cw  configuration. 
Indeed,  if  the  energy  is  delivered  to  the  target  in  a sequence  of  short 
powerful  bursts,  the  heat  energy  deposited  in  the  atmosphere  by  the  previous 
pulses  may  have  time  to  be  carried  away  from  the  propagation  path  by  conduction, 
free  convection  or  forced  convection  before  the  next  pulse  is  emitted. 

Therefore,  favourable  conditions  of  application  may  be  obtained  through  a 
judicious  choice  of  the  laser  parameters.  A parametric  study  of  the 
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computer  solutions  is  performed  for  a collimated  Gaussian  laser  beam 

in  the  presence  of  a constant  cross  wind.  It  shows  that  not  only 

can  the  adverse  effects  of  thermal  blooming  be  eliminated  but  that  an  actual 

enhancement  of  the  average  power  density  on  the  target , with  respect  to  the 

non-blooming  situation,  can  be  achieved. 

Section  2.0  presents  a brief  review  of  the  theoretical  background 
and  of  the  working  hypotheses.  It  also  discusses  the  model  of  multipulse 
heating  and  defines  the  similarity  parameters.  The  numerical  method  is 
described  in  Section  3.0  and  Section  4.0  gives  and  analyzes  the  results 
for  a collimated  Gaussian  beam.  The  work  was  performed  at  DREV  between 
January  and  June  1976  under  PCN  34C04  (formerly  PCN  07E04,  Project  No. 
97-01-38),  Propagation  of  Laser  Beams. 


2-0  THEORETICAL  BACKGROUND 

The  propagation  model  used  to  construct  our  algorithm  is 
fully  described  in  Ref.  7.  It  differs  from  more  conventional  mo- 
dels in  that  the  wave  equation  is  broken  into  two  equations:  one  equation 

for  the  geometrical  rays  and  one  equation  for  the  complex  amplitude  of  the 
field  defined  on  the  geometrical  phase  front.  This  separation  is  exact  and 
has  the  practical  advantage  of  allowing  longer  propagation  steps  in  the 
numerical  solution  of  the  finite  difference  version  of  the  governing 
equations . 

2.1  Working  Hypotheses  and  Governing  Equations 

Our  study  examines  the  thermal  blooming  phenomenon  under 
dominant  forced  convection  cooling,  i.e.  when  the  heat  energy  supplied  by 
absorption  is  totally  balanced  by  forced  convection  losses.  The  velocity 
component  of  the  wind  (and/or  of  the  beam  motion)  in  the  direction  transverse 
to  the  optical  axis  is  assumed  constant.  Only  linear  changes  in  the  thermo- 
dynamic properties  of  the  medium  are  considered  and  compressibility  effects 


UNCLASSIFIED 

3 


and  kinetic  cooling  are  neglected.  Finally,  longitudinal  derivatives  are 
neglected  compared  with  crosswise  derivatives.  Under  these  conditions,  the 
medium  equations  are  reduced  to  a single  temperature  equation  which  has  the 
f orm: 

39  W30  a T, 

3t  + 3?  = — I(X’y’Z;t)’  (1) 

o p 

where  9 is  the  temperature  fluctuation,  W the  wind  speed,  a the  absorption 
coefficient  of  the  medium  at  the  laser  beam  wavelength,  pq  the  ambient  density 
and  Cp  the  specific  heat  at  constant  pressure  of  the  medium  and  I(x,y,z;t) 
is  the  beam  irradiance. 


We  consider  a laser  source  pulsating  at  a single  and 
constant  angular  frequency  m,  neglect  polarization  and  make  the  paraxial 
approximation.  If  the  scalar  electric  field  is  written  in  the  form: 


E = Aexp[jk(z+4>) -jaz-jwt]  , 
the  wave  equation  becomes,  after  separation: 

(lr  * v-vt)v  - ■ 


(L-^v.vJa  - -S*VY  ♦ k'' 


V 2a, 


where  j=/-l , V = 7^ , 7^  is  the  gradient  operator  in  the  transverse 
plane,  n and  nQ  respectively  the  perturbed  and  unperturbed  index  of 

refraction,  and  k the  wave  number  defined  as.  follows:  k = n w/c. 

o 

Equation  (3)  is  the  eikonal  equation  and  thus  V represents  the  vector 
angle  subtended  by  the  geometrical  rays.  Equation  (4)  is  the  equation 
for  the  complex  amplitude  A defined  on  the  geometrical  phase  front: 


tor  tne  compi 

(fr  ‘ v-7<) 


indicates  derivation  along  a geometrical  ray,  jAV^.V 


represents  the  geometrical  effects  of  the  converging  and  diverging 

rays  and  j 7 2A  is  the  diffractional  contribution. 

2k  z 


FIGURE  la  - 


Schematic  illustration  of  the  time  evolution  of  the 
irradiance  function  I(x,y,z;t)  of  a multipulse  laser 
beam. 
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FIGURE  lb 


Schematic  illustration  of  the  steady  state  time  evolu- 
tion of  the  periodic  temperature  distribution  that 
results  from  the  heating  of  the  air  by  a multipulse 
laser  beam. 
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2 . 2 Steady  State  Multipulse  Heating 

Multipulse  thermal  blooming  has  been  studied  at  other  laboratories 
e.g.  Refs.  2,  3,  4,  5,  11,  12  and  13.  However, the  solution  we  propose 
for  the  temperature  profile  is  written  in  a somewhat  more  convenient 
form  that  does  not  involve  summation.  Therefore  it  is  worth  reviewing 
its  derivation  in  some  details. 

We  are  searching  for  a steady  state  solution  to  the  tem- 
perature equation  (1)  when  the  forcing  function  I is  a periodic 
sequence  of  short  pulses.  The  function  I(x,y,z;t)  is  schematically 
illustrated  in  Fig.  la.  The  pulse  duration  b is  assumed  negligible 
compared  with  the  period  of  pulsation  T,  i.e.: 

b « T.  (5) 

The  resulting  steady  state  temperature  solution  is  a periodic  function 
of  period  T;  it  is  schematically  represented  in  Fig.  lb. 

For  the  present  purpose,  we  restrict  ourselves  to  a propa- 
gation regime  in  which  single-pulse  thermal  blooming  can  be  neglected. 

This  approximation  is  valid  if  the  pulse  duration  is  smaller  than  the 
thermalization  time  of  the  molecular  energy  levels  excited  by  absorption, 
or  smaller  than  the  transit  time  of  sound  waves  across  the  beam.  Typically, 
these  time  scales  are  of  the  order  of  1 ms.  Moreover,  the  incom- 
pressibility hypothesis  of  Section  2.1  requires  that  the  period  of  pul- 
sation T be  much  longer  than  these  characteristic  time  scales,  which  is 
consistent  with  relation  (5).  Physically,  these  approximations  mean 
that  the  properties  of  the  medium  seen  by  the  m^  pulse  are  completely 
determined  by  the  incompressible  heating  that  results  from  the  absorption 
of  the  (m-1)  pulses  which  have  preceded  it.  Therefore,  since  b <<  T, 
it  is  clear  from  Figs,  la  and  lb  that  we  need  only  to  solve  the  temper- 
ature level  attained  just  prior  to  the  sharp  increase  that  follows 
the  arrival  of  the  mth  pulse,  where  m is  large  enough  to  have  steady 


state. 
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Between  pulses,  equation  (1)  for  the  temperature  profile  takes 


the  form: 


39  W39 

rr-  + tt—  =0. 

3t  3y 


A general  solution  to  equation  (6)  is  given  by: 

0[x,y,z;t]  = 0[x,y-W(t-tQ) ,z] . (7) 

Integrating  equation  (1)  between  t * mT-b/2  and  t = mT+b/2,  we  obtain 
under  the  approximation  expressed  by  relation  (5) : 

6[x,y,z;mT+b/2]  - 9[x,y,z;mT-b/2]  = D[x,y,z]  , (8) 

o P 


where 


mT+b/2 


D[x,y,z]  = /l[x,y,z;t]dt 

•'.rr*  V 


mT-b/2 


is  the  energy  density  profile  of  each  pulse. 

Since  the  overall  solution  is  periodic,  we  have  from 
equation  (7) : 

0[x,y,z;mT+x]  = 0 [x ,y-Wx , z ;mT+b/2] , (10) 

where  b <_  t <_  T-b  and  where  0[x,y,z;mT+b/2]  is  the  temperature 
distribution  immediately  following  the  firing  of  the  mt^1  pulse  which 
serves  as  initial  condition.  In  particular,  we  have  from  equation  (10) 
at  x = T-b/2 : 


0[x,y,z;(m+l)T-b/2]  = 6[x,y-WT,z;mT+b/2] , (11) 

where  Wb/2  in  the  spatial  coordinate  of  the  right-hand  side  function 
of  equation  (11)  was  neglected  in  accordance  with  relation  (5). 
Replacing  (m+l)T  by  mT  in  the  left-hand  side  function  of  equation  (11), 
which  is  justified  since  the  solution  for  9 is  periodic,  and  using 
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equation  (8)  to  evaluate  the  temperature  jump  between  (mT-b/2) 
and  (mT+b/2)  in  the  right-hand  side  of  equation  (11),  we  obtain  the 
following  solution  for  the  steady  state  temperature  distribution  just 
prior  to  the  arrival  of  the  in  pulse: 

9[x,y,z;mT-b/2]  = 0 [x,y-WT,z;mT-b/2] 

+ rV  D[x,y-WT,z],  (12) 

o p 

and  the  boundary  condition  is  simply: 

lim  0[x,y,z;mT-b/2]  = 0.  (13) 

y-^-oo 

The  (-00)-  limit  is  actually  the  upwind  boundary  of  the  laser  beam. 

In  summary,  equation  (12)  together  with  boundary  condition 
(13)  constitutes  a simple  recurrence  formula  to  compute  the  temperature 
distribution  seen  by  the  m^  pulse.  In  the  approximation  of  no  single- 
pulse thermal  blooming,  this  is  sufficient  to  solve  the  steady  state 
average  irradiance  distribution  of  repetitively  pulsed,  high  average 
power  laser  beams.  Note  that  solution  (12)  contains  no  summation 
over  the  (m-1)  earlier  pulses  in  contrast  with  expressions  given  in 
Refs.  3,  11  and  12. 

2 . 3 Scaling  Laws  and  Similarity  Parameters 

For  numerical  calculations  to  be  practical,  it  is  important 
to  define  and  use  the  scaling  laws  and  the  similarity  parameters  of  the 
problem.  This  minimizes  the  number  of  parameters  that  have  to  be 
varied  to  obtain  a full  picture  of  the  situation.  The  technique 
essentially  consists  in  writing  the  governing  equations  in  dimension- 
less form  and  the  non-dimensional  parameters  emerging  from  the  process 
constitute  the  similarity  parameters.  This  analysis  for  the  present 
propagation  model  is  fully  described  in  Refs.  7-8  and,  therefore, 
the  results  will  be  simply  recalled  here. 
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The  time  average  irradiance  profile  of  the  multipulse  laser 
beam  is  defined  as  follows: 


<I>  = AA*e 

t 


-az 


(14) 


where  the  amplitude  A satisfies  equation  (4).  If  P designates  the  time 
average  power  and  a, some  representative  radius  of  the  beam,  the  char- 
acteristic scale  for  the  wave  amplitude  A is  given  by: 


Ac  = (P/ira2)2, 

the  characteristic  scale  for  the  temperature  variation  6 by: 
aPT 


9 


p C na2 
o p 


the  characteristic  scale  for  the  angle  of  the  geometrical  rays 
V by: 


f ootP  M 

rpocPn°aW)  ’ 


the  longitudinal  length  scale  by: 

1 

Up  C n a3W 
1 o p ~ 


!ip  C n a3W  )2 

-°  P-° \ 

aaP  | 


and  the  crosswise  length  scale  by: 


t = a. 
c 


(15) 


(16) 


(17) 


(18) 


(19) 


n 


Using  the  scaling  laws  defined  by  equations  (15-19)  and 
denoting  the  non-dimensional  functions  by  a superscript  ' , one  obtains 
the  following  governing  equations: 

e'U,e,n;mT-b/2]  = e - [c,5-2<,n;mT-b/2] 

+ (A'A'*)  [C,C-2K,n;mT-b/2] .exp(-Yn) , (20) 


(|_  .V.v;)y. 


= 7;e' 


(21) 
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(k*v-7;)A'  ■ -»A'7{-?' 

j V2A’ 
+ *—  t 

2F  1 

> 

(22) 

where  the  non-dimensional  variables  are 

defined  as 

follows : 

n = z/£p, 

(23) 

£ = x/l  , and 
c 

(24) 

^ = y/V 

(25) 

From  examination  of  the  non-dimensional  governing  equations 
(20-22),  we  have  the  following  three  similarity  parameters: 

WT 

K = 2iT  ’ 

(26) 

Y = , and 

P 

(27) 

F = ka2 

(28) 

r~ 

p 


For  a collimated  beam  of  given  initial  geometry,  k, 
y and  F are  the  only  numbers  that  are  required  to  completely 
specify  the  problem  of  steady  state  multipulse  thermal  blooming 
under  the  conditions  and  approximations  listed  in  the  preceding 
paragraphs,  k is  called  the  displacement  number  which  gives  an 
estimate  of  the  number  of  beam  diameters  a heated  particle  of 
air  is  transversely  drifted  between  successive  pulses,  y is  the 
absorption  number  related  to  the  effect  of  overall  linear  de- 
pletion and  finally  F is  the  Fresnel  number  which  measures  the 
magnitude  of  cw  thermal  blooming  distortion  with  respect  to 
natural  divergence  or  diffractional  effects. 
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3.0  NUMERICAL  METHOD 


The  non-dimensional  governing  equations  (20-22)  are 
integrated  numerically.  The  algorithm  and  computational  details 
are  fully  described  in  Ref.  10.  The  system  of  finite  difference 
equations  is  implicit,  and  the  nonlinearity  in  equation  (21) 
and  the  coupling  between  equations  (20-22)  are  solved  through 
iterations.  The  predictor-corrector  approach  of  Ref.  10  has  been 
replaced  by  a simpler  and  unconditionally  stable  method  of  alter- 
nating directions  (Ref.  14) . 


The  functions  0' 


V'  and  A'  are  symmetric  with  restect  to 


the  5-  axis  while  is  antisymmetric.  All  functions  are  defined  on 
a 30  x 60  homogeneous  grid  and  the  mesh  size  is  adjustable  at  each 
propagation  step  to  optimize  the  utilization  of  the  grid  points. 

The  longitudinal  increment  An  is  recalculated  at  each  step  to 
satisfy  the  accuracy  criterion.  Typically,  An  is  of  the  order 
0.1,  which  corresponds  in  practical  applications  to  100-300  m.  These 
relatively  long  forward  computational  steps  are  possible  since  the 
wave  amplitude  A is  defined  on  the  geometrical  phase  front  which,  at 
high  Fresnel  numbers,  constitutes  almost  the  entire  thermal  blooming 
phase  distortion. 


4.0  RESULTS  FOR  COLLIMATED  GAUSSIAN  LASER  BEAMS 


Sample  results  of  the  predicted  irradiance  distribution  of  an 
originally  Gaussian  laser  beam  with  F = 20  and  y = 0.45  are  illustrated 
in  Figs.  2-5.  Figs.  2-3  give  the  continuous  wave  (cw)  thermal 
blooming  irradiance  profile  (k  = 0)  at  n = 1.15  in  the  form  of  iso-irra- 
diance  contours  and  isometric  projection  respectively.  Similar  results 
at  n = 1.12  for  multipulse  thermal  blooming  (k  = 0.8)  are  reproduced  in 

Figs.  4-5. 
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Figs.  2-5  show  that  there  are  significant  differences  between 
cw  and  multipulse  thermal  blooming.  The  time  average  peak  irradiance  of 
the  multipulse  distribution  is  nearly  twice  as  large  as  that  of  the  cw 
distribution  and  approximately  1.9  times  the  maximum  irradiance  of  the  non- 
bloomed  beam.  The  point  of  maximum  irradiance  of  the  multipulse  pattern 
is  located  almost  on  the  line  or  sight  whereas  it  is  displaced  upwind  by 
0.9  beam  radius  in  the  cw  case.  The  beam  divergence  in  the  direction 
normal  to  the  wind  is  near  the  diffraction  limit  in  the  multipulse  regime 
but  exceeds  this  limit  in  the  cw  regime.  Physically,  the  absorbed  heat 
in  the  multipulse  situation  is  shifted  on  the  downwind  side  of  the  beam 
when  the  m^  pulse  is  emitted.  Therefore,  the  upwind  half  of  the  beam 
sees  no  refractive  index  gradient  and  propagates  unperturbed  whereas  the 
downwind  half  sees  a gradient  which  tends  to  bend  its  rays  into  the  wind. 
Hence,  conditions  may  be  achieved  where  the  two  portions  of  the  beam  are 
almost  exactly  superposed  which  results  in  a factor  of  2 enhancement  in  the 
power  density  delivered  to  a target.  The  results  depicted  in  Figs.  4-5 
approach  this  ideal  situation:  the  iso-irradiance  contours  show  that 

a large  fraction  of  the  available  laser  power  is  concentrated  in  the  up- 
wind half  of  what  would  be  an  un-bloomed  beam  distribution. 

The  results  presented  in  Figs.  2-5  would  approximately 
correspond  to  the  irradiance  patterns,  at  3 km,  of  a 1 MW  (time  average) 
collimated  Gaussian  laser  beam  of  e-folding  diameter  (irradiance 
profile)  equal  to  60  cm.  The  environmental  conditions  would  be:  sea  level, 

absorption  coefficient  equal  to  0.17  km  ^ and  a lateral  wind  velocity  equal 
to  10  m/s. 


In  summary,  Figs.  2-5  show  that  multipulse  thermal  blooming  may 
create  favourable  conditions  of  application  leading  to  an  intensification 
of  power  density  on  the  target  which  are  worth  investigating. 


-t.l 
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FIGURE  2 - Iso-irradiance  contour  lines  of  the  computed  beam  irradiance 
pattern  at  n = 1.15  for  F = 20,  y = 0.45  and  k = 0 (cw) . 

The  irradiance  levels  are  defined  by  steps  of  0.1  between 
0.1  and  0.9  times  the  peak  irradiance. 


FIGURE  3 


Computed  beam  irradiance  profile  at  n = 1.15  for  F 

v = 0.45  and  k = 0 (cw) . 


M FIGURE  4 - Iso-irradiance  contour  lines  of  the  computed  beam  irradiance 

pattern  at  n = 1.12  for  F = 20,  y = 0.45  and  k = 0.8 
. > [multipulse) . See  Fig.  2 for  a definition  of  the  irradiance 

K--  levels. 


FIGURE  5 


Computed  beam  irradiance  profile  at  n = 1.12  for  F 
Y = 0.45  and  tc  = 0.8  (multipulse). 
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4.1  Experimental  Verification 


Very  few  experimental  data  on  multipulse  thermal  blooming  are  yet 
available.  The  only  direct  measurements  known  to  the  author  were  per- 
formed in  a simulation  experiment  conducted  by  Buser  and  co-workers 
(Ref.  13).  A train  of  pulses  from  a TEA-C02  laser  (typically  0.5  J/pulse) 
was  propagated  through  a 1-m  long  laterally  translated  absorption  cell. 

A low  power  cw-C02  probe  beam  was  used  coaxially  to  monitor  the  dis- 
tortion of  the  irradiance  pattern.  Of  special  interest  to  us  are  the 
results  illustrated  in  Fig.  9 of  Ref.  13.  This  figure  shows  the  evolution 
of  the  probe  beam  peak  irradiance,  relative  to  the  non-blooming  case, 
as  a function  of  time  following  the  first  pulse  of  the  collinear 
repetitively  pulsed  laser  beam.  The  data  run  for  the  first  three  pulses. 
Given  the  conditions  of  application,  especially  k = 0.69,  steady  state 
should  be  approximately  attained  after  two  pulses.  The  measured  relative 
peak  irradiance  at  the  time  of  firing  of  the  third  pulse  is  somewhere 
between  1.30  and  1.35. 

Exact  numerical  modelling  of  the  experiment  described  in  the 
preceding  paragraph  is  not  possible  with  the  present  version  of  the 
program.  This  is  so  because  the  absorbing  medium  in  that  experiment  is 
pure  C02  at  a pressure  of  7 atm  which  saturates  at  the  high  peak  power  of 
the  multipulse  beam  but  does  not  at  the  low  power  of  the  cw  probe  beam  which 
is  actually  used  to  make  the  measurements.  The  unsaturated  absorption 
coefficient  is  0.42  m”1  and  the  extinction  coefficient  of  the  multipulse 
TEA  laser  is  0.22  m”  . Therefore,  the  differences  are  significant.  However, 
we  can  account  for  this  effect  in  the  following  approximate  manner.  In  the 
calculation  of  the  scaling  length  and  of  the  similarity  parameter  F, 
we  use  the  unsaturated  value  of  the  absorption  coefficient  since,  at 
< = 0.69,  the  contribution  to  thermal  blooming  deformation  comes  mainly 
from  absorption  on  the  wings  of  the  beam  where  saturation  does  not  occur 
or  is  less  pronounced.  But,  to  evaluate  the  absorption  number  y,  we 
use  the  saturated  coefficient  since  this  number  is  related  to  linear 
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depletion  from  which  the  value  of  0.22  m * was 
obtain  for  the  conditions  of  Fig.  9 of  Ref.  13 
for  the  similarity  parameters: 

tc  = 0.69  , 

Y “ 0.35  , 

F = 3.91  , 

and  the  data  are  reported  for  n = 0.63.  If  these  numbers  are  used  in  the 
computer  program,  we  find  for  the  relative  steady  state  peak  irradiance 
of  the  probe  beam  at  n = 0.63,  just  prior  to  the  firing  of  a pulse  from 
the  TEA  laser,  a value  equal  to  1.35.  Therefore,  the  agreement  with  the 
reported  data  is  almost  exact.  This  excellent  agreement  is  very  gratifying 
when  compared  with  results  of  other  workers.  For  example,  it  is  found  in 
Ref.  13  that  the  relative  peak  irradiance  intensification  predicted  by 
the  numerical  code  of  Ref.  11  is  about  50%  smaller  than  the  experimental 
one. 

4 . 2 Parametric  Study 

In  view  of  the  predicted  multipulse  thermal  blooming  enhancement 
shown  in  Figs.  2-5  which  is  confirmed  experimentally  by  the  data  of 
Ref.  13,  it  is  important  to  examine  the  dependence  of  the  resulting 
intensification  upon  the  similarity  parameters  k,  F and  y.  This  constitutes 
the  most  practical  way  of  defining  optimal  operating  conditions. 

A simple  measure  of  beam  quality  is  provided  by  the  re- 
lative time  average  peak  irradiance.  The  adjective  relative  is  used  in 
this  document  to  designate  normalization  of  a variable  with  respect 
to  the  same  variable  calculated  from  the  linearly  absorbed  but 
un -bloomed  irradiance  pattern.  The  relative,  steady  state,  time 
average  peak  irradiance  is  plotted  in  Fig.  6 as  a function  of  the 
displacement  number  k for  a fixed  Fresnel  number  F = 20  and  a 


determined.  We  thus 
the  following  values 
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fixed  absorption  number  y = 0.45.  Four  curves  are  presented 
corresponding  to  four  non-dimensional  propagation  distances  n = 0.4, 

0.8,  1.0  and  1.1  respectively.  All  curves  are  shown  to  pass  through 
a maximum  point  for  a value  of  < located  somewhere  between  0.7  and 
0.9.  Fig.  7 illustrates  the  same  results  for  a different  choice  of 
F and  y,  namely  F = 50  and  y=  0.2.  Qualitatively,  the  curves  shown 
in  both  figures  resemble  one  another  very  closely  except  that  the 
enhancement  is  generally  higher  for  the  conditions  of  Fig.  7.  Optimal 
peak  irradiance  intensification  at  n = 1.1  is  as  large  as  95%  for 
F = 50,  y = 0.2  and  85%  for  F = 20,  y = 0.45.  The  optimum  displacement 
number  in  all  instances  is  approximately  < = 0.8.  Therefore,  it  can 

be  stated  as  a first  conclusion  that  positive  or  favourable  steady- 
state  thermal  lensing  can  be  expected  in  almost  all  operating  conditions 
if  the  repetition  rate  is  chosen  to  give  < - 0.8. 

Since  k = 0.8  represents  optimal  intensification,  it  is 
interesting  to  study  how  the  latter  varies  with  respect  to  F and  y. 

The  relative  peak  irradiance  at  ranges  of  n = 0.4,  0.8,  1.0  and  1.1 
is  plotted  in  Fig.  8 versus  Fresnel  number  for  < = 0.8  and  y = 0.2. 
Intensification  steadily  increases  with  F until  F = 20,  where  it  passes 
through  a maximum  point  at  the  longer  propagation  distances,  and  then 
settles  to  what  appears  to  be  a nearly  constant  value  for  F > 30. 
Therefore,  at  propagation  ranges  near  n = 1,  Fig.  8 reveals  the  exis- 
tence of  an  optimal  Fresnel  number  whose  value  is  slightly  greater  than 
20;  for  example,  the  curve  n = 1.1  shows  a peak  irradiance  inten- 
sification greater  than  100%  for  20  < F < 30.  Finally,  there  remains  to 
examine  how  the  intensification  under  the  most  favourable  conditions 
defined  by  k = 0.8  and  F = 20  changes  with  the  third  and  last  similarity 
parameter  y.  This  is  done  in  Fig.  9 which  shows  that  the  optimal 
gain  slowly  decreases  with  increasing  y for  n = 1.0  and  1.1  but  remains 
nearly  constant  for  n = 0.4  and  0.8.  The  largest  drop  in  intensification 
is  only  20%  over  y ranging  from  0.1  to  1.0.  Hence,  it  appears  that  the 
relative  optimal  peak  irradiance  is  only  slightly  affected  by  linear 
absorption.  This  result  is  very  interesting  since  the  absorption 
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number  y is  the  one  parameter  that  cannot  be  controlled  once  the  operating 
wavelength  or  wavelengths  have  been  chosen. 


Although  the  time  average  peak  value  of  the  transverse  irradiance 

distribution  constitutes  a reasonably  useful  criterion  to  appraise  the  beam 

quality,  it  gives  insufficient  information  since  the  thermally  bloomed  beam 

has  a very  complex  and  non-self-similar  geometry.  However,  this  deficiency 

can  be  avoided  by  measuring  the  beam  quality  in  terms  of  the  time-space  average 

irradiance,  denoted  <I>  , over  the  area  of  the  transverse  plane  which 

’ t-s’  v 

intercepts  a given  fraction  of  the  total  available  power.  To  compute  this 
spatial  average,  we  need  to  calculate  the  area  delimited  by  the  iso-irra- 
diance  lines  of  the  irradiance  pattern  as  those  illustrated  in  Figs.  2 and 
4.  A computer  program  was  written  that  routinely  performs  this  calculation 
for  an  arbitrary  number  of  iso-irradiance  contour  lines  of  arbitrary  shape. 

The  spatial  averaging  inside  the  N*"*1  level  contour  line  is  simply  done  as 
follows : 


<I>. 


t-s 


2S‘ 


1 T 
+ <J>t 


] [S1  - Si_1], 


(29) 


where  S1  designates  the  area  bounded  by  the  i*  contour  curve  and  <I>^  is 

0 r 

the  corresponding  irradiance  level;  <I>  is  the  time  average  peak  irradiance 
0 X 

and  S = 0.  The  chosen  irradiance  levels  are  consecutively  ordered  from  the 

highest  to  the  lowest  value.  If  M denotes  the  rank  of  the  zero-irradiance 
level,  the  fraction  f(N)  of  the  total  power  which  passes  inside  the  N1"*1- 
level  contour  line  is  given  by: 


f(N) 


<I> 


N 

t-s 


<I> 


M 

t-s 


(30) 


The  results  for  the  time-space  average  irradiance  are  presented 

in  Figs.  10-13.  The  parametric  analysis  is  performed  in  the  same  manner 
as  in  Figs.  6-9.  Two  series  of  curves  are  given,  distinguished  by  a suffix 
a or  b following  the  figure  number.  Series  a and  b refer  respectively  to  the 
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relative  time-space  average  irradiance  over  the  area  receiving  63%  and  85% 

63% 

of  total  power.  These  quantities  are  denoted  as  follows:  DCI <I>° 

KfcL  t - S 

DC.<I>°  °,  63%  representing  the  fraction  of  total  power  comprised  within 
the  circle  of  e-folding  radius  of  a Gaussian  beam  and  85%  being,  to  the 
nearest  5%,  the  fraction  of  total  power  contained  inside  the  circle 
bounded  by  the  first  dark  ring  of  the  diffraction  pattern  of  a uniformly 
illuminated  circular  aperture.  These  quantities  are  truly  independent 
of  beam  shape  and  therefore  they  should  constitute  very  practical  and 
useful  criteria  of  laser  performance  under  thermal  blooming  distortion. 

63% 

Qualitatively,  the  results  shown  in  Figs.  10-13  for  _ .<I> 

85%  KuL.  t — s 

and  DCI<I>°  ° are  very  similar  to  those  given  in  Figs.  6-9  for  the 

relative  time  average  peak  irradiance.  First,  there  is  a definite 

optimal  displacement  number  k and  its  value  obtained  from  Figs.  10  and 

11  is  also  near  k = 0. 8-1.0  as  it  was  found  from  the  peak  irradiance  data. 

The  thermal  blooming  enhancement  of  the  time-space  average  irradiance, 

although  smaller  than  for  the  peak  irradiance,  is  nevertheless  substantial. 

For  example,  at  n = 1.1,  F = 50  and  y = 0.2,  it  is  better  than  50%  for 

nt,.<I>  and  better  than  40%  for  nci<I>°  °.  Secondly,  the  curves  of 

KhL  t "S  Kt L L ~ S 

constant  n in  Figs.  12a  and  12b  show  that  the  time-space  average  irradiance 
gain  monotonically  increases  with  Fresnel  number  F before  reaching  a 
constant  asymptotic  value.  Therefore,  there  is  no  clear  maximum  as 
previously  observed  in  Fig.  8 with  regard  to  the  relative  peak  irradiance, 
but  one  can  still  define  a favourable  Fresnel  number  range.  Indeed, 
examination  of  Figs.  13a  and  13b  reveals  that  F > 20  ensures  optimal 
thermal  blooming  gain  for  propagation  distances  n ± 1.1.  Finally, 

Figs.  13a  and  13b  show,  as  before,  that  the  absorption  number  y has  only 
a mild  influence  on  multipulse  thermal  blooming  intensification. 
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FIGURE  6 - Relative  maximum  values  of  the  irradiance  profiles  at  distances 

n = 0.4,  0.8,  1.0  and  1.1  as  functions  of  similarity  parameter 
k for  fixed  F = 20  and  y = 0.45.  Relative  refers  to  the  irra- 
diance distribution  of  the  un-bloomed  but  linearly  absorbed 
beam. 
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FIGURE  8 - Relative  maximum  values  of  the  irradiance  profiles  at 

distances  n = 0.4,  0.8,  1.0  and  1.1  as  functions  of 
similarity  parameter  F for  fixed  ic  = 0.8  and  y = 0.2. 
See  Fig.  6 for  a definition  of  the  term  relative. 


FIGURE  9 - Near  optimal  (<  = 0.8,  F = 20)  relative  maximum  values 

of  the  irradiance  profiles  at  distances  n = 0.4,  0.8, 
1.0  and  1.1  as  functions  of  similarity  parameter  y. 

See  Fig.  6 for  a definition  of  the  term  relative. 
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5.0  CONCLUSIONS 


Numerical  calculations  of  steady  state  irradiance  distributions  • 

of  repetitively  pulsed,  high  average  power,  collimated  Uaussian  laser  beams 
show  that  multipulse  thermal  blooming  can  constitute  a positive  lens. 

Intensifications  of  the  time  average  peak  irradiance  better  than  100%  and 
of  the  time-space  average  irradiance  (over  t*''*  area  containing  85%  of  total 
power)  better  than  40%  were  demonstrated  numerically.  These  numbers  are  at 
least  twice  as  high  as  those  quoted  by  other  workers  (Refs.  5 and  13)  but 
available  experimental  data  (Ref.  13)  confirm  almost  exactly  our  predictions 
whereas  those  of  Ref.  11  (as  reported  in  Fig.  9 of  Ref.  13)  are  too  low  by 
approximately  the  same  factor  of  two. 

A dimensional  analysis  of  the  governing  equations  has  yielded 
three  similarity  parameters:  1)  the  displacement  number  x which  is  the  * 

reciprocal  of  the  number  of  pulses  per  flow  time,  2)  the  Fresnel  number 
F which  measures  the  magnitude  of  cw  thermal  blooming  distortion  with 
respect  to  diffractional  divergence,  and  3)  the  number  y related  to  linear 
power  losses  by  absorption.  These  three  numbe-s  completely  specify  the 
problem  of  steady  state  multipulse  thermal  blooming  of  collimated  beams. 

A parametric  study  of  computer  solutions  for  collimated  Gaussian  beams 
has  demonstrated  that  optimal  intensified  power  density  on  the  target  can 
be  obtained  under  conditions  defined  by: 

0 . 8 < x < 1 . 0 , and 

F > 20 

The  remaining  parameter  y has  no  optimal  range  but  it  should  be  chosen  as 
small  as  possible  both  to  minimize  linear  absorption  and  to  increase  thermal 
blooming  intensification.  However,  its  influence  on  the  latter  phenomenon  (F, 
x and  n being  the  same)  is  rather  mild. 


These  conclusions  reveal  a favourable  or  advantageous  aspect  of 
thermal  blooming  in  the  form  of  multipulse  rather  than  cw  operation  which 
certainly  must  be  included  in  the  design  of  high-power  military  laser  systems. 
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Of  course,  this  report  examines  only  one  aspect  of  the  problem; 
abstraction  was  made  of  technical  feasibility  and  of  other  linear  and  non- 
linear effects.  To  mention  only  the  most  important,  these  effects  are 
particle  scattering,  platform  jitter,  atmospheric  turbulence,  single-pulse 
thermal  blooming  and  gas  breakdown.  Some  of  these  effects  are  in  competition. 
For  example,  pulse  duration  should  be  chosen  short  enough  to  avoid  adverse 
single-pulse  thermal  blooming  but,  at  the  same  time,  long  enough  to  prevent 
gas  breakdown.  Ref.  5 studies  the  interrelation  of  all  these  propagation 
phenomena  in  the  light  of  a simplified  empirical  model,  but  additional 
work  is  still  required  to  account  for  all  effects  in  a more  general  and  precise 
manner.  Our  immediate  plans  at  DREV  are  to  incorporate  the  scintillation 
effects  of  turbulence  in  the  computer  algorithm. 
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intensite  et  le  nombre  y relie  J I'absorption  lin£aire.  Une  etude  parametrique  en  intensity  et  le  nombre  y rel i£  Z I'absorption  lin^aire.  line  ^tude  parametrique  en 

fonction  de  ces  trois  param&tres  demontre  que  des  conditions  optimales  de  gain  fonction  de  ces  trois  paramfctres  d£montre  que  des  conditions  optimales  de  gain 

dans  le  cas  d'un  faisceau  gaussien  collimate  existent  pour:  0.8  s x s 1.0  et  F j 20.  dans  le  cas  d'un  faisceau  gaussien  collimat€  existent  pour:  0.8  s r s 1.0  et  F ; 

L 'ampl  if  icat  ion  decroft  generalement  avec  l'accroissement  de  y mais  de  fa«,on  peu  L 'ampl  ification  d^croft  generalement  avec  1 'accroissement  de  y mais  de  faqon  peu 

sensible.  sensible. 


